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In this work it is investigated if fermionic sources could be responsible for accelerated periods 
during the evolution of a universe where a matter field would answer for the decelerated period. The 
self-interaction potential of the fermionic field is considered as a function of the scalar and pseudo- 
scalar invariants. Irreversible processes of energy transfer between the matter and gravitational 
fields are also considered. It is shown that the fermionic field could behave like an infiaton field in 
the early universe and as dark energy for an old universe. 
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I. INTRODUCTION 

The search for constituents responsible for accelerated 
periods in the evolution of the universe is a fundamen- 
tal topic in cosmology. Usually the formulations include 
elements of general relativity, field theory and thermo- 
dynamics, putting under analysis the evolution of space- 
time variables like the scale factor, its acceleration and 
energy densities. Several candidates has been tested for 
describing both the inflationary period and the present 
accelerated era: scalar fields, exotic equations of state 
and cosmological constants. 

Another possibility is to consider fermionic fields 
as gravitational sources for an expanding universe. 
Fermionic sources has been investigated using sev- 
eral approaches, with results including exact solutions, 
anisotropy-to-isotropy scenarios and cyclic cosmologies 
(see, for example P,|3,I3)- 

In the present work the connection between general rel- 
ativity and the Dirac equation is done via the tetrad for- 
malism, where the components of the tetrad or "vierbein" 
play the role of the gravitational degrees of freedom. 
The interactions between the constituents are modeled 
through the presence of a non-equilibrium pressure term 
in the source's energy-momentum tensor. Besides, it is 
considered a self-interaction term for the fermionic con- 
stituent, in the form of a potential that can assume sev- 
eral forms (included the Nambu-Jona-Lasinio case 3). 

While testing fermionic sources as responsible of ac- 
celerated periods different regimes are possible. In a 
young universe scenario the fermion produces a fast ex- 
pansion where matter (included via a barotropic equation 
of state) is created till it starts to predominate and the 
initial accelerated expansion gives place to a decelerated 
era, dominated by a matter field, which ends when the 
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fermionic field predominates again leading to an acceler- 
ated era. In this case the fermionic field plays the role of 
the infiaton in the early period of the universe and of dark 
energy for the old universe, without the need of a cosmo- 
logical constant term or a scalar field. In an old universe 
scenario an initially matter dominated period gradually 
turns into a dark (fermionic) energy period when an ac- 
celerated regime starts and remains for the rest of the 
evolution of the system. 

The manuscript is structured as follows. In section 
II we make a brief review of the elements of the tetrad 
formalism used to include fermionic and matter fields in 
a dynamical curved space-time. The field equations for 
an isotropic, homogeneous and spatially flat universe are 
derived in section III . In section IV it is presented the 
analysis of the different scenarios in which the fermionic 
constituent answers for accelerated eras and the transi- 
tions from accelerated to decelerated periods and vice- 
versa. Finally, in section V we display our conclusions. 
The metric signature used is (-)-,—,—,—) and units have 
been chosen so that SttG = c = h= k=l. 



II. DIRAC AND EINSTEIN EQUATIONS 

In this section we present briefly the techniques that 
are used to include fermionic sources in the Einstein the- 
ory of gravitation and for a more detailed analysis the 
reader is referred to |E IS S 111- The starting point is 
that the gauge group of general relativity does not ad- 
mit a spinor representation and the tretad formalism is 
invoked to solve the problem. Following the general co- 
variance principle, a connection between the tetrad and 
the metric tensor g^i, is established through the relation 

9^.l ^ e^e^^Tjab, a = 0,1, 2, 3 (1) 

where ej^ denotes the tetrad or "vierbein" and rjab is the 
Minkowski metric tensor. Here Latin indices refer to the 
local inertial frame whereas Greek indices to the general 
system. 
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As it was said above, the main objective of this work 
is to describe the behavior of fermions in the presence of 
a gravitational field, and the next step is to construct an 
action for this system. The Dirac lagrangian density in 
Minkowski space-time is 

Ld='^ bPl^dai' - ida^h^i^] - - V, (2) 

where the spinors are treated as classically commuting 
fields Above, m is the fermionic mass, ?/)=-(/!t--yO ^^g, 
notes the adjoint spinor field and the term V , which is 
a function of -0 and i/;, describes the potential density of 
self-interaction between fermions. The general covariance 
principle imposes that the Dirac-Pauli matrices must 
be replaced by their generahzed counterparts = 6^7°, 
whereas the generalized Dirac-Pauli matrices satisfy the 
extended Clifford algebra, i.e., {r'',r''} = 2^'^". 

In a second step we need to substitute the ordinary 
derivatives by their covariant versions 

(3) 

where the spin connection fl^ is given by 

n^ = -\9,AKx-^UdAWj\ (4) 

with V^^ denoting the Christoffcl symbols. Hence, the 
generally covariant Dirac lagrangian becomes 

Ld = ^[^T^'D^^ - {D^4>)V^^\ - - V. (5) 

The field equations are obtained from the total action 

S{g,i',il') = J V^Ltd^x, (6) 

where Lt = Lg + Ljj + Lm is the total lagrangian den- 
sity. Lg — R/2, with R denoting the curvature scalar is 
the gravitational Einstein lagrangian density for fermions 
which are minimally coupled to the gravitational field. 
Lu is the Dirac lagrangian density (jsj and is the 
lagrangian density of the matter field. 

From the lagrangian density (O, through Euler- 
Lagrange equations, we obtain the Dirac equations for 
the spinor field and its adjoint coupled to the gravita- 
tional field 

iVD^i^-m^-^ = 0, iD^i^T^'+mlp+— = 0. (7) 
dip dip 

The variation of the action ^ with respect to the 
tetrad leads to Einstein field equations 

R^l.v — -^gfj,uR = —Tfj,^, (8) 

where Tp,y is the total energy-momentum tensor which is 
a sum of the energy-momentum tensor of the fermionic 
field T'^" and of the matter field T^^^" , i.e., Tt"" = T^'' + 



Tj^^. Furthermore, the symmetric form of the energy- 
momentum tensor of the fermionic field is given by 

T'f" = - r^F'^D''^/; + v-r^iJ^V' 

/ 4 

-D^i^ViP - D'^lpViP] - g^'^Lo. (9) 

III. FIELD EQUATIONS 

The Robertson- Walker metric incorporates the homo- 
geneity and isotropy hypotheses for the universe. Here we 
consider a spatially flat universe described by the metric 

ds^ = dt^ ~ a{tf{dx^ + dy^ + dz'^), (10) 

where a{t) refers to the cosmic scale factor. 

The total energy-momentum tensor for an isotropic 
and homogeneous universe - which is composed by 
fermionic and matter fields and where dissipative pro- 
cesses are taken into account - is written as 

(T^^) = diag(p, -p -w,-p- w, ~p - w). (11) 

Above, the total energy density p and the total pres- 
sure p are given as a sum of the corresponding terms 
of the fermionic and matter fields, i.e., p = Pf + Pm 
and p = Pf + Pm- Moreover, the quantity za refers to 
a non-equilibrium pressure which is related to dissipative 
processes during the evolution of the universe and repre- 
sents an irreversible process of energy transfer between 
the matter and the gravitational field 9] . 

Thanks to the Bianchi identities, the covariant dif- 
ferentiation of Einstein field equations l|HJl leads to the 
conservation law of the total energy-momentum tensor 
T^'^-i, = 0, hence it follows by using the representation 
(|11|1 the evolution equation for the total energy density: 

p + 3H{p + p + w) ^0, (12) 

where the dot refers to a differentiation with respect to 
time and H = d{t)/a{t) denotes the Hubble parameter. 
Furthermore, from Einstein field equations l(HJl it follows 
the Friedman and acceleration equations 

H^ = \p. - = -i(p + 3p + 3n7), (13) 
3 a b 

respectively. Only two equations from (|12|l and H13|) are 
linearly independent. 

For the metric H1U|I the tetrad components read 

eS = «, er = ^*.". (14) 

Also, the Dirac matrices become 

F° - 7", F* = 4t7\ = -i^T^T^T^T^ = 7^ 
a{t) 

(15) 
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from which the spin connection components are obtained, 
yielding 



^0 = 0, = ^d{t)^V- 



(16) 



For an isotropic and homogeneous universe the 
fermionic field is an exclusive function of time. So the 
Dirac equations iQ become 

iP + ^HiP + im'y°iJj + i'y°^ = 0, (17) 
2 dtp 



^ + -HliJ - im^-i° - i—-f° = 0, (18) 
2 dip 

thanks to equations (|15|l and (|16(l . 

The non-vanishing components of the energy- 
momentum tensor of the fermionic field follow from Q 
together with lO, - yielding 



(19) 



(T.)\^(r,)V(T,)V.-^|-||, (20) 

which are only functions of "0 and Tp. By identifying 
the components of the energy-momentum tensor of the 
fermionic field as (Tf)'^^ = diag(p/, — p/, — p/, —pf), one 
can obtain from l|17|) and H18|l by the use of 119|) and 
(|20|l the conservation law for the energy density of the 
fermionic field: 



Pf + 3H{pf+pf) = 0. 



(21) 



Hence, the evolution equation for the energy density of 
the fermionic field decouples from the energy density of 
the matter field, and we have from (|12|l and (|21|) : 



(22) 



where the term —ZHw could be interpreted as the energy 
density production rate of the matter field (see e.g. 

According to the extended (causal or second-order) 
thermodynamic theory the non-equilibrium pressure 
w{t) obeys an evolution equation, whose linearized form 
reads 



3r,H, 



(23) 



where r denotes a characteristic time and r/ is the so- 
called bulk viscosity coefhcient. 



IV. COSMOLOGICAL SOLUTIONS 

In order to analyze cosmological solutions for the 
field equations of the previous section, we have to spec- 
ify the potential density of self-interaction between the 



fermions V. According to the Pauli-Fierz theorem V 
is an exclusive function of the scalar invariant (ipip)'^ 
and on the pseudo-scalar invariant (itpj^ip)'^, i.e., V — 
V{{ipip)'^, (tipj^ip)'^). Here we are interested in analyzing 
self-interaction potentials of the form 



F = A[/3i(i^V)' + /32(«V'7V)' 



(24) 



where the coupling constant A is a non-negative quantity 
and n is a constant exponent. We shall analyze three 
cases, namely, (i) (3i = 1 and 02 = where is a function 
only of the scalar invariant; (ii) Pi — and /?2 = 1 where 
V depends only on the pseudo-scalar invariant and (iii) 
Pi — 02 — ^ where ^ is a combination of the scalar 
and pseudo-scalar invariants. The Nambu-Jona-Lasinio 
potential is represented by the last case with n = 1. 

The energy density and the pressure of the fermions 
for the potential (|24|l are given by 

Pf = m(^V') + A [Pii^tPf + /32(«V^7V)']" , (25) 



Pf = {2n-l)X [Pi (V-V)" + 02 " , (26) 

respectively, thanks to (|20l) . 

We infer from (|26|l that the fermions could be classi- 
fied according to the value of the exponent n. Indeed, for 
n > 1/2 the fermions represent a matter field with posi- 
tive pressure (n > 1/2) or a pressureless fluid (n = 1/2), 
whereas for n < 1/2 the pressure of the fermions is neg- 
ative and they could represent either the inflaton or the 
dark energy. 

For massless fermions, the pressure is connected with 
the energy density by a simple barotropic equation of 
state p/ — (2n—l)pf and it follows from the conservation 
equation (|21|) for the fermions that pj cx 1/a^". We 
shall not analyze this case here, since the behavior of 
the fermionic field does not differ from that of a matter 
field when n > 1 /2 or from that of a bosonic field when 
n < 1/2. Moreover, for the massive case, we shall deal 
with only the case where the fermionic field behaves as 
infiaton or dark energy, i.e., the case where n < 1/2. 

For the pressure of the matter field we shall adopt a 
barotropic equation of state, i.e., Pm — WmPm 

with < 

Wm 1- Furthermore, the coefficient of bulk viscosity 77 
and the characteristic time r are assumed to be related 
with the energy density p hy rj — ap and r = ry/p, where 
a is a constant Q. 

The system of field equations we shall investigate in 
order to find the cosmological solutions are: 
(a) the acceleration equation 



a _ 1 
a 6 



(p/ + Pm + 3p/ + 3p„ + 3tn); (27) 



(b) the evolution equation for the energy density of the 
matter field 



3i/(p„ 



Pr, 



0; 



(28) 
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(c) the evolution equation for the non-equihbrium pres- 
sure 

TW + m ^ -i-qH; (29) 

(d) the Dirac equation 117|) , which in terms of the spinor 
components ip = (V'l; V'2, "03, "04)"^, becomes 




- 21(0-1 V'l +i^\i^2 - 01^3 - 0204) 




V* = 0. (30) 



In equation (|15|l we have introduced the following abbre- 
viations 



(9(V'0)2 



dV 



(31) 



9(075-0)2 

Equations (|27|l through (|30() consist of a system of 
seven coupled ordinary differential equations for the fields 
a{t), pm{t), T^{t), 01 (i), 02 (i), tp3{t) and ip4{t) and in the 
following subsections we shall find solutions of this sys- 
tem of equations for given initial conditions. 



A. Accelerated-decelerated-accelerated regime 

Let us first analyze the case that corresponds to the 
evolution of the early universe, where the fermionic field 
plays at the beginning the role of an inflaton and the 
matter is created by an irreversible process through the 
presence of a non-equilibrium pressure w. The initial 
conditions we have chosen for t — (by adjusting clocks) 
are 



a(0) = l, 0i(O) = O.lz, 02(0) = 1, 

MO)^i, Pm(0)-0, tI7(0)=0. 



03(0) = 0.3, 
. (32) 

The last two initial conditions correspond to a vanishing 
energy density of the matter field and a vanishing non- 
equilibrium pressure at t = 0. The conditions chosen here 
characterize qualitatively an initial proportion between 
the constituents in the corresponding era (i.e., in this 
case we have a predominating fermionic field over the 
matter, indicating the initial inflationary state) 

Since equation H27() is a second-order differential equa- 
tion we need to specify an initial condition for d(0). This 
condition follows from the Friedmann equation (|13|l i . i.e.. 



a(0) = a(0) 



P/(0)+p™(0) 




(33) 



FIG. 1: Acceleration field a vs. time t. 



with pf{0) determined from equation H25|) . However, 
we have to specify some parameters in order to obtain 
numerical solutions of the coupled system of differential 
equations (|27|) through H30I) . These parameters are: (a) 
^1 Pii P2 and n which define the self-interacting potential 
(I24|l : (b) m which is related to the mass of the fermionic 
field; (c) Wm which defines the matter field through its 
barotropic equation of state pm = WmPm and (d) a which 
is connected with the bulk viscosity term rj = ap. In 
order to plot the figures 1 and 2 we have chosen the fol- 
lowing values for these parameters: 

A = 0.1, /3i=/32 = l, n = 0.3, 
771 = 0.01, w„i = 1/3, a = 1.0, and a — 1.2, 

(34) 

which correspond to a fermionic field with a negative 
pressure, described by a self-interacting potential that 
depends on the scalar and pseudo-scalar invariants and 
a matter field of massless particles that could describe a 
radiation field. 

In figure 1 it is plotted the acceleration field d as func- 
tion of time t for two different values of a = 1.0 and 1.2, 
whereas in figure 2 it is shown the behavior of the energy 
densities of the fermionic p f and matter p,„ fields as func- 
tions of time t. We infer from these figures that there ex- 
ists an accelerated period where the fermionic field domi- 
nates and the matter field is created at the expenses of an 
irreversible process of energy transfer between the mat- 
ter and the gravitational field. The accelerated period 
is followed by a decelerated era which is dominated by 
the matter field. Due to fact that the self-interaction po- 
tential tends to a constant value for large values of time, 
the energy density of the fermionic field overcomes the 
energy density of the matter field and the universe goes 
into another accelerated period. It is noteworthy that 
for the bosonic case where the potentials are exponen- 
tials or inverse power-laws, one has to add a constant 
value ~ which is similar to introduce a cosmological con- 
stant term - in order to return to an accelerated era after 
the accelerated-decelerated period (see, for example [lO|). 
Here the same self-potential interaction which plays the 
role of an inflaton field at the beginning plays the role 
of a cosmological constant for later times and could be 
identify as dark energy. We note from the figures that 
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FIG. 2: Energy densities of fermionic pf and matter pm fields 
vs. time t. 




FIG. 3: Acceleration field a vs. time t. Dashed line - a = O.f ; 
straight line - a = 0.05; dotted line - without irreversible 
processes. 



for large values of the coefficient a it follows that: (a) 
the energy density of fermionic field decays more rapidly 
causing a larger accelerated period and (b) the energy 
density of the matter field has a more significant growth 
and leads to a larger decelerated period. 

The same conclusions above could be obtained for 
a self-interacting potential which is only a function of 
the scalar invariant (/3i = 1, /32 = 0) or of the scalar 
pseudo-invariant (/3i — (32 — 1). However, all cases 
are strongly dependent on the exponent n of the self- 
interacting potential and one can obtain different behav- 
iors in which there exits only an accelerated period for the 
universe or the universe begins with a decelerated period. 
This last case will be analyzed in the next subsection. 

It is worth mentioning that the behavior of the fields 
found in our analysis is not restricted to the initial con- 
ditions given in ^62\ and the values for the parameters 
given in l|34() . In fact, these represent typical values that 
describe qualitatively the transitions under investigation 
in this case, i.e., the accelerated-decelerated case when 
matter emerges as the predominating constituent and a 
second situation where dark energy, represented by the 
fermionic field, overcomes the matter field and leads to a 
final accelerated period. 



B. Decelerated-accelerated regime 

Let us now investigate an old universe dominated by 
non-relativistic matter or dust {wm = so that pm — 0) 
where a fermionic field is present but with an energy den- 
sity smaller than that of the matter field. We shall con- 
sider the same initial conditions (by adjusting clocks) as 
those of previous subsection for a(0) and V'i(O) through 
■04(0), but for energy density of the matter field we regard 
it as twice that of the fermionic field, i.e., /O„i(0) — 2/9/(0). 
With respect to the parameters, we shall choose the same 
values as above for X, (3i, P2, 'm and n, but different values 
for a, namely a = 0.1 and 0.05. We have also considered 
the case where the irreversible processes during the evo- 
lution of the universe are absent. For this last case the 
evolution equation of the non-equilibrium pressure H29|) 
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FIG. 4: Energy densities of fermionic pf and matter pm fields 
vs. time t. Dotted and dotted-dashed lines - without irre- 
versible processes; straight and dashed lines - a = 0.1. 



was not taken into account. The acceleration and the 
energy densities as functions of time are plotted in fig- 
ures 3 and 4, respectively. We observe from these figures 
that there exists a transition from a high deceleration - 
where the universe is dominated by the non-relativistic 
matter field - to a small acceleration - where the uni- 
verse is dominated by the fermionic field which plays the 
role of dark energy. By considering irreversible processes 
there is no sensitive change in the energy density of the 
fermionic field, but the matter field decays more slowly 
so that the accelerated period begins later. 

The same conclusions of last subsection regarding the 
self-interaction potential are also valid for this case. 



V. FINAL REMARKS AND CONCLUSIONS 

One question that could be formulated is about the 
importance of the irreversible processes during the evo- 
lution of the universe. To answer this question one should 
discuss the relevance of the bulk viscosity coefficient in 
the different periods of the universe's evolution. In the 
above calculations it has been used a dimensionless coef- 
ficient, but if one goes back to its dimensional expression, 
it is found that it has to be multiplied by a factor pro- 
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portional to the initial value of the Hubble parameter. 
Since wc know that the present value of that parameter 
is very small, this coefficient would be playing an impor- 
tant role only for early periods of the universe. On the 
other hand for the late period, although its contribution 
is very small, the presence of the viscosity is necessary 
for describing the thermodynamic dissipative effects in 
an expanding universe. 

We have investigated the possibility that a fermionic 
field - with a self-interacting potential that depends on 



the scalar and pseudo-scalar invariants - could be the 
responsible for accelerated regimes in the evolution of 
the universe. We have shown that the fermionic field 
behaves like an inflaton field for the early universe and 
later on as a dark energy field, whereas the matter field 
was created by an irreversible process connected with a 
non-equilibrium pressure. Moreover, for an old deceler- 
ated universe dominated by non-relativistic matter the 
fermionic field plays again the role of dark energy and 
drives the universe to an accelerated regime. 
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